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Abstract 

We analyze theoretically the influence of backward propagating waves on the primary threshold 
when a linearly polarized light impinges at normal incidence on a planarly aligned nematic liquid 
crystal films. We show, that the primary threshold, as a function of the phase delay induced by 
the nematic layer, exhibits oscillations. The amplitude of oscillations depends strongly on the drop 
of the refr activity indices of the nematic and outer media at the boundaries. 

PACS numbers: 42.70.Df, 05.45.-a, 42.65. Sf 
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I. INTRODUCTION 



Many aspects of complex light-induced nonlinear orientational phenomena in liquid crys- 
tals have been investigated during the last few decades, with a particular emphasis on the 



mc 



associated with the elas- 



]. From the optical point 



dynamical richness arising from the strong light-matter coup 
tic and anisotropic optical properties of ordered mesophases 
of view, nematics are uniaxial media with the optical axis along the local average direc- 
tion of the molecular axis called director, n. Depending on the boundary conditions and 
bulk ordering characteristics, a light-driven orientational instability may take place above a 
threshold intensity usually referred to as the optical Freedericksz transition (OFT). Besides 
their fundamental interest, the optically induced phenomena in liquid crystals seem to be 
finding their way for technological applications. It has been recently proposed 5K7) to real- 
ize all-optical photonic switching devices using dielectric periodic structures infiltrated with 
homeotropically or planar aligned nematics. 

When considering a linearly polarized light field impinging at normal incidence onto a 
planar nematic cell, one should distinguish between two situations when the polarization 
plane is parallel or perpendicular to the initial director uniform distribution. In both cases, 
twisted elastic reorientation modes appear above a light intensity threshold. The first one 
corresponds to the coherent polarization conversion process mediated by two-beam coupling 
between the incident extraordinary wave and a lower-frequency noisy ordinary component. 
The latter field component is due to scattering of the incident light by the director orienta- 
lional fluctuations and experiences gain during the propagation of light throughout the cell 
s), 1^]. The second situation, which corresponds to an ordinary polarized incident light, is 
associated with the twist Freedericksz transition. Up to now, this particular case has only 
been considered theoretically because of the predicted threshold for orientational instability 



is turned out to be very high 12] . Recently, the twist OFT has been revisited and it was 
shown that the stationary bifurcation solution derived in Ref. 12] is valid for optically thin 
nematic layers. Indeed, for thicker films, the twisted OFT takes place via a Hopf bifurcation 
[lo| . which has not been demonstrated experimentally yet. It is worth noting that the planar 
geometry has attracted an interest following seemingly experimentally accessible intensity 



threshold using a periodic structure 



111- 



In this work, in an attempt to guide future experimental investigations we quantitatively 
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estimate the influence of unavoidable back and forward light scattering on the twist OFT 
threshold that are due to boundary refractive indices mismatch in any real sample. This 
effect is usually neglected, which is well justified for homeotropic alignment. However, when 
planar alignment is considered, the perfectly matching condition can be achieved only by 
utilizing anisotropic substrates. This is obviously an inconvenient solution in practice. Our 
analysis reveals the crucial role of backward waves in the calculation of the threshold itself. 



II. THE MODEL 



We consider a linearly-polarized plane wave with wavelength A impinging at normal 
incidence (along z) on a planarly aligned nematic liquid crystal layer of a thickness L (see 
Fig. [p. The unperturbed director n lies along the x axis and the light is polarized along 
the y direction, i.e. we deal with an incident ordinary wave. Theoretical description of 
optically induced orientational phenomena in nematics is a well-posed problem that consists 
of coupled nematodynamic and Maxwell's equations for the light propagation inside the 
anisotropic medium |4| . There are several common simplifications that are employed to solve 
these equations. First of all, fluid flow associated with the director reorientation (backflow 



effect) is neglected. Indeed, flow p 
slight quantitative differences only 



ays only a passive role and in most known cases lead to 

In 

131 . H4J . In this study we neglect all transverse effects and 
assume that all variables depend solely on z. It allows to simplify drastically the Maxwell's 
equations. Therefore, the light beam might be treated as a plane wave and any experiment 
aiming to verify the theoretical predictions should be performed by using a laser beam with 
a diameter much larger than the thickness of the nematic layer. The applicability of such a 
one-dimensional approximation has already been validated for homeotropic nematic samples, 
|z. It turned out that it also works rather well for the beam diameters of the order 



with 


n 


z. I 


of L 


15, 


16| 



16]. Importantly, the realization of spatially inhomogeneous transverse structures 



in the (x, y) plane is hardly probable. As was shown in 17] in the limit of thin layers, it 
is a case only for strongly anisotropic ratio between elastic constants which is not a case 
for common nematics. Next, we assume that the director always keeps planar orientation 
at the boundaries (strong anchoring conditions), i.e. n(z = 0,t) = n(z = L,t) = n . 
We introduce the representation adapted to our geometry in terms of the spherical angles 
Q(z,t) and t) such that n = (cos $ sin O, sin $ sin O, cos O). However, as follows from 
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the simulations, the director reorients in the (x, z) plane only for the all the regimes explored 
in the present study. Therefore O = 7r/2 and the director is described by the twist angle 
$ only. Such a situation is intuitively expected because the light should transmit a lot of 
energy to deform the director across the layer. The strong boundary conditions in terms of 
the twist angle read 

$(z = 0,*) = $(* = 7T,t) = 0. (1) 



III. PRIMARY INSTABILITY WITH IDEAL LIGHT FIELD BOUNDARY CON- 
DITIONS 



We first neglect the effect of reflection at the boundaries z = 0, 7r (ideally matched 
conditions) and calculate the thresholds by performing the linear stability analysis of the 
planar state. We eventually obtain the following integro-differential equation for $ [lO ] 

<9 f $ = <9 z 2 $ + 2pA 2 ( $ + A ^ $(z') sin[A(z' - z))dz'\ (2) 

where p = I/I c is the dimensionless incident light intensity with I c = 8n 2 cK2n e Sn/[X 2 (n e + 
n )] and 

ttA = 2nL5n/\ (3) 

is the phase delay between e- and owaves through the whole layer with birefringence 5n = 
n e — n a , where n e and n are refractive indices of the o and e waves, respectively. Here we 
introduce the normalized length z — > ztt/L and time t — > t/r where r = 71L 2 / '{t^K?) is 
a characteristic relaxation time with the rotational viscosity 71 and the twist Frank elastic 

n 

constant K 2 |4|. 

We perform then the linear stability analysis writing $(2, t) = 3>(z) exp(crt) and inserting 
the expansion 

$0, t) = <Pm(t) sin(mz). (4) 

m 

into Eq. (j2J) with further projection on cp„ (Galerkin method). The resulting eigenvalue 



problem is solved then numerically (see [lCj] for details). 

The analysis shows that the primary instability of the planar state is a stationary bi- 
furcation only below a critical value A c = 0.64 of the phase delay and a Hopf bifurcation 



otherwise, as indicated by the blue dashed line in Fig. |2j We notice that the symmetry S is 
spontaneously broken at the threshold in former case. The chan ge f rom stationary to Hopf 
bifurcation is accompanied by a jump of the threshold intensity [10 ]. 



IV. PRIMARY INSTABILITY WITH REALISTIC LIGHT FIELD BOUNDARY 
CONDITIONS 

As a next step, we consider realistic boundary conditions at z — 0, 7T, namely an 
isotropic dielectric material that sandwiches the nematic layer with refractive index n on f 
For convenience, we introduce the labels i and t for the first (z < 0, incident side) and 
the second (z > n, transmitted side) semi-infinite surrounding dielectric material, and 
iV for the nematic layer. The modified linear equation requires a derivation of the ex- 
pressions for the o and e waves inside nematic to the first order in $ that account for 
backward waves contributions. This is done by using the Berremen's formalism [18 ], 
when the Maxwell equations are cast in the matrix form for the four component vector 
1 F T = (E x , H y , E y , —H x ). Here E(z, t), H(z, t) are the amplitudes of the electric and mag- 
netic fields, E(r,t) = l/2(E(«,*)e _iwt + c.c), H(r,t) = l/2(H(z, t)e~ iuJt + c.c), that vary 
slowly in time compared to u) . The vector \l/ then can be written as a superposition of the 
four proper electromagnetic waves (two forward- and two backward-propagating) propagat- 
ing across the layer (k x = k y = 0) with the same value u for i, t and iV media (Oldano's 



formalism 



201]). namely 
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* = E h¥ = M l + + + M 4 = T0, (5) 

i=i 

where T is the matrix whose columns are eigenvectors ^ and <fi is the vector of amplitudes 
f = (A, /a, / 8 , / 4 )- 

The vector of amplitudes <fi at the left side of the first boundary plane (z = 0") and at 
the right side of the second boundary plane (z = tt + ) can be written as 

T (O") = (0, oq, r e , To), t (tt+) = (t e , t , 0, 0) (6) 

The scattering problem consists of expressing the amplitudes r e , r D , t e , t D of the reflected 
and transmitted e- and o waves as a function of the amplitude a of the incident o-wave. 
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The continuity of the tangential components at z = 0,n of the vectors E and H requires 



^((T) = ^(0 + ), *(tt - ) = *(?r + ). 

so that we can write 

T^(O-) = T JV (O + )0(O + ),T w (7r-)0(7r-) = T t( f>(n+), 
where Tj and T t are as follows 



(7) 



(8) 
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(9) 



Due to the boundary conditions for the director we have T N (0 + ) = T N (n ) = Tft given by 

/— — n \ 



T° - — 

N ~ V2 



/n e 
fn~ e 




-y/n^ 
1 
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(10) 



\ v /n^ -y/n J 

The evolution of the electromagnetic field inside nematic can be written by introducing the 
propagation matrix K(n) as 



0(tt-) = K(tt)0(O + ), 



(11) 



where 



K(tt) = iafc e i7rfcon ° x 
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(12) 



Here 



a = 



n 2 e -n 2 



_ n e + n 
^ 5n ' 



(13) 



and 

Ci = / e- iAz $(z)dz, f 2 = / e- iAqz $(z)dz. (14) 
Jo Jo 

Finally, a set of four independent equations that allows to express r e , r Q , t e and t as a 

function of ao is derived 

0(0-) = Tr^^K-^Tr)^)- 1 ^^). (15) 

In general, the resulting expressions for the reflected and transmitted waves r e ,r Q ,t e and t 
are cumbersome, however, it is instructive to consider the particular case n ont = n Q that 
simplifies the algebra a lot and gives 

ia qA £ [e*^*-*) + qe ~ iA( i^~ z )] $(z) dz 

r Q = 0, (16) 
za gAe2 wA (^i) [ ge - jA2 + e iAqz ] $(z) dz 



t 



g2g— iirqA gj-7rA 



where g = (n e + n )/5n. As is expected, r D = because an incident ordinary wave is 
perfectly matched (n out = n ). Thus the o wave is entirely transmitted acquiring only the 
phase according to Eq. f[T6"j) . 

Next, we determine the light field distribution inside nematic by integrating the equations 
for the amplitudes fi(z), since the initial condition 0(O + ) = (T^)- 1 T i (0, a , r e , r ) T is 
found. Finally, the linearized expression for the electromagnetic term can then be calculated, 
leading to the following integro-differential equation for $ 

<9 t $ = d 2 z <$> + 2pA 2 f $ + A jf sin[A(z' - z)]dz' + A • Im[F(z)]J , (17) 

where F(z) describes the contribution of the backward waves. Then, the stability analysis 
procedure is similar to the case of ideal boundary conditions (see Sec. IIIII) . In the particular 
case n out = n Q , the latter function writes explicitly 

rz ( p iAz I -iqAz\ rn 

F{z) = / e-^~ z ^(z')dz' + {6 mA + ge > ■ \ [e^^ + q e~ iA ^} *(z) dz (18) 
Jo e ~ Q e Jo 

The results are summarized in Fig. [2] where primary instability threshold is plotted in 

the (A,p) plane with (solid curves, for the two cases n out = n Q and n out = 1) and without 



(dashed curve) taking into account the contribution of backward waves generated at the 
nematic boundaries. The presence of backward waves has a clear qualitative effect, namely, 
the oscillation of the threshold with respect to the phase delay A of the nematic layer. Such 
an oscillation is driven by the dependence of the function F(z) on A. In the case n ont = n , 
the resulting oscillation has a period 2/q [see Eq. (ITS]) ]. It corresponds to a period of 
A/(n e + n Q ) when considering the dependence of the threshold versus thickness L. The fact 
of oscillations itself emphasizes the interference between forward and backward propagating 
waves to be at the origin of the modulation of the threshold. Accordingly, the modulation 
depth drastically increases when the refractive index contrast at the nematic boundaries is 
increased, as shown when taking n out = 1 (see Fig. [2]). To conclude this section, let us note 
that in a practical situation the representative case will correspond to n out = 1. Indeed, 
a sample is typically made of a nematic slab sandwiched into two glass substrates (with 
refractive index n glass ~ n G ) surrounded with air, hence we expect the air/glass interfaces to 
control the modulation of the OFT threshold, which therefore corresponds qualitatively to 
the case n ont = 1. 

V. CONCLUSIONS 

We report on the effect of back- and forward scattering on the Freedericksz transition 
in a nematic layer with planar alignment. The light intensity and the phase delay induced 
by the layer are the control parameters of the system. We demonstrate that the primary 
threshold exhibits oscillations with the phase delay. This effect is especially important when 
the jump of refractivity index at the boundaries of nematic layer occurs. 
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FIG. 1: Light-matter interaction geometry: linearly polarized light along the y direction incident 
perpendicularly on a nematic layer with the director no || x (planar state). The components of the 
director n are described in terms of the twist angle 
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FIG. 2: (Color online) Primary instability threshold of the planar state on the (A, p) plane with 
(solid lines) and without (blue dashed line) the effect of backward waves. Orange (light gray) line: 
refractive index of outer media n ou t = n Q . Red (dark gray) line: n out = 1. Lines with A < A c 
(A > A c ) correspond to stationary (Hopf) bifurcation. 
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